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ABSTRACT 



Relativistic Cowling approximation in which all metric perturbations are omitted is applied to 
non-axisymmctric infinitesimal oscillations of uniformly rotating general relativistic polytropes. 

Frequencies of lower order f-modes which are important in analysis of secular instability 
driven by gravitational radiation are investigated and neutral points of the mode along 
equilibrium sequences of rotating polytropes are determined. Since this approximation becomes 
more accurate as stars are more relativistic and/or as they rotates more rapidly, we will be able 
to analyze how a rotation period of a neutron star may be limited by this instability. 

Possible errors in determining neutral points caused by omitting metric perturbations are 
also estimated. 



Subject headings: relativity — stars:rotation — stars:oscillation 



Introduction 



Since the discovery of non-axisymmetric secular instability of rotating stars driven by gravitational 
radiation (Chandrasekhar 1970, Friedman fc Schutz 1978j , Friedman 1978), its possibility to limit a 



rotational period of a neutron star has been pointed out by several authors. As the instability sets in at 
neutral points in the absence of viscosity, it would become the main issue to determine these neutral points 
for rapidly rotating stars in general relativity. However, thus far, only several investigations in Newtonian 
gravity ([Managan 1985 , Imamura, Friedman & Duriscn 1985) and one post-Newtonian investigation 
( putler fc Lindblom 199% have been made, because it has been very difficult to manage general relativistic 
perturbations of rotating stars completely. 

In such a situation, Yoshida & Kojima (1997) have examined the accuracy of the relativistic version 
of the Cowling approximation for oscillations of slowly rotating relativistic stars. Comparing the results 



by the relativistic Cowling approximation (hereafter RCA; for earlier literatures, see e.g. McDermott et 



al. 1983, Finn 1988, Lindblom & Splinter 1990, Ipser & Lindblom 1992) with 'exact' eigenfrequencies 
obtained by fully relativistic analysis, they argued that the approximation works well for slowly rotating 
relativistic stars. This suggests that we can determine with a fairly good accuracy the neutral points by 
this approximation. 
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We here investigate the "counter-rotating" f-modc oscillations of rapidly rotating relativistic stellar 
models by the RCA. Though, as seen in the Newtonian results ( Robe 1968| ), errors of f-mode eigenfrequencies 
are larger than those of p- or g- modes, absolute values of errors arising from the RCA are likely to become 
smaller for relativistic and rotating stars compared to the Newtonian values (see below). Furthermore we 
will estimate roughly the amount of the corrections to the results of the RCA. 



2. Basic Equations and Numerical Method 

Here we investigate only uniformly rotating polytropic stars. Polytropic relation is defined by 
(p/e c ) — K(e/e c ) 1+1 / JV , where e c is the central energy density, K is a measure of strength of gravity of the 
star and other symbols have their usual meanings. 

The stationary axisymmctric equilibrium states are obtained by the KEH scheme (Komatsu, Eriguchi, 
& Hachisu 1989). The space coordinates (r,0,(p) and the time coordinate t are used in this paper. As a 
parameter characterizing the stellar rotation we adopt the ratio of the rotational energy to the gravitational 



energy T/\W\ (see Komatsu et al. 1989) 



In the present study, all the Eulerian perturbations of the metric components are omitted. Thus 
our basic equations consist of (1) the perturbation equation of the baryon number conservation and (2) 
the perturbation equation of the energy momentum conservation. The perturbed Eulerian variables are 
assumed to have harmonic dependences on the time and the azimuthal angular coordinate: ~ e -wi+im<^ 
where m is an integer and a is the frequency. The boundary condition of the problem is that the Lagrangian 
perturbation of the pressure vanishes at the original stellar surface. Since the basic equations closely 
resemble those of Newtonian gravity, they are solved by the same scheme developed for the Newtonian 
perturbations (Yoshida fc Eriguchi 1995). R 



3. Results 

3.1. Test Calculations 

We have checked our RCA code by applying it to two cases, i.e. spherical stars in general relativity 
and rotating stars in Newtonian gravity. First our frequencies for relativistic spherical stars agree well with 
those of other RCA calculations of Yoshida & Kojima (1997). As for rapidly rotating Newtonian stars, 
by comparing our results with those of the Newtonian Cowling approximation which are obtained by the 
modified code of Yoshida & Eriguchi (1995), we have found that our RCA code can reproduce the neutral 
points of f-modes with enough accuracy, i.e. within errors of 2%. 



2 Although Ipser & Lindblom (1992) have presented an elegant scalar potential formalism, we do not use that formalism here 
because we prefer a formulation which is parallel to that for the Newtonian case (Yoshida and Eriguchi 1995). We thus derive 
the relativistic version of the linearized continuity equation and of three components of the linearized Euler equation and solve 
them directly. See appendix for explicit form of the equations. 
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3.2. Neutral Points along Sequences 

From intersection of T/\ W\ — a curves of equilibrium sequences and <r = 0we can obtain neutral points 
which are tabulated in Table 1. 

In the framework of Newtonian gravity critical values of T 1 /! V4^| for the m = 2 'bar' mode have been 
considered to be universal, i.e. r/lH^I ~ 0.14 irrespective of compressibility and/or the rotation laws 
( Tassoul 197§| ). However, even for N = 0.5 polytropes we could not find a neutral point against this 
bar mode in weak gravity limit, i.e. for small n. This may results from the crudeness of the Cowling 
approximation for lower order modes. In the Newtonian Cowling approximation, smaller frequencies are 
obtained than those calculated from the full treatment. As a consequence, the critical values of T/|VF| at 
neutral points tend to be raised in this approximation. Thus even if N < 0.808 neutral points for the bar 
mode may not be reached before the termination of equilibrium sequence due to mass-shedding. 

On the other hand, neutral points for the bar mode could be found for stronger gravity, although the 
value of T/|Vy| does not coincide with the Newtonian universal value 0.14. Remarkable is that even for 
soft equations of state (N > 0.808) with sufficiently strong gravity, neutral points for the bar mode do 
appear. This is not, however, so surprising because general relativistic effect on equilibrium states becomes 
significant. [|] 

Our numerical computational results can be summarized as follows: (a) for sequences with N and k 
kept fixed, values of T/|VF| at neutral points become smaller as the value of m increases; (b) for sequences 
with N and m kept fixed, values of T/|W| at neutral points become smaller as the value of K increases; (c) 
for sequences with m and K kept fixed, values of T/\ W\ at neutral points becomes smaller as the value of N 
increases. 

In Newtonian gravity, normalized eigenfrequency scales as: ^^ vp ~ ^~p~ > wnere Pc is the mass density 
at the centre, p is the mean density of the star and I is the zenithal quantum number. From this relation, 
normalized frequencies become smaller for models with higher mass concentration toward the centre. 
Because of strong gravity, matter distributions are concentrated towards the central region for relativistic 
models. In other words, relativistic stars can be considered to become softer effectively. Consequently 
frequencies become lower for larger values of n. 



4. Discussion How Far is the RCA Reliable ? 

Before applying the RCA to realistic neutron stars, it is important to note how far this approximation 
can be reliably applied. In particular we need to estimate how the strength of relativity and the amount of 
rotation affect the results. 

It should be noted that eigenfunctions of energy density perturbation for f-modes are peaked near 
the surface. It implies that surface regions are important for f-modes. On the other hand, as is known 
in Newtonian models, relative importance of gravitational perturbations reduces as the amount of mass 
participating in oscillations decreases. 

Consequently relativistic gravity can be considered to improve the approximation because the effective 



3 After submitting this paper, we received a preprint of Stergioulas & Friedman (1997), where they also report to have 
obtained neutral points of m = 2 mode for softer equations of state with taking metric perturbations into consideration. 
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concentration of the matter due to relativity makes the role of metric perturbations less important. 
Furthermore perturbed density distribution of rotating models is more sharply peaked toward the surface 
of the star as the stars are spun up. Therefore the accuracy of the Cowling approximation can be relatively 
increased for rotating models than spherical ones. Thus these two factors can be said qualitatively to have 
advantages in deciding neutral points. 

However, since there have been no investigations on perturbations of rotating stars in full general 
relativity, the influence of rotation on the errors due to the RCA cannot be quantitatively estimated by 
direct comparison with the results of full calculations. Thus we need some device to have quantitative 
estimation about the errors. 

Our proposal is as follows. We assume that qualitative dependence of the mode on stellar rotation in 
general relativity is similar to that of the Newtonian counterpart. In Newtonian theory, differences between 
the eigenfrequencies obtained by the Cowling approximation and those by the full analysis including 
gravitational perturbations decrease monotonically as stars are spun up. If this tendency is the same for 
general relativistic rotating stars, the maximal error of the RCA in determining neutral points can be 
estimated. 

Let us select one T/|W| — a curve for a certain equilibrium sequence computed by the RCA. Since the 
maximum error of this curve is assumed to be that of the spherical model, the amount of the maximum 
error can be calculated by comparing our results with exact values of spherical models in general relativity. 
A parallel displacement of this T/|W| — a curve by this amount of the maximum error yields a curve with 
the maximal error in determination of neutral points (Fig. 1). Since the amount of displacement required 
to obtain the correct eigenfrequency from the value by the Cowling approximation decreases as the value 
of T/|W| is increased, the 'true eigenfrequency curve' with gravitational perturbations passes between the 
curves of the Cowling approximation and of the curve of the maximal error defined above. 

Therefore if the error of the eigenfrequency for the spherical model and the inclination of T/|VF| — a 
curve of the RCA at the approximated neutral point are known, the maximal error in determining a neutral 
point will be determined by: 



where A<7o is the maximum error of the frequency. 

For N — 1 polytropes we have obtained the maximal errors of T/|W| at neutral points for several 
modes (Fig. 2). We also display the post-Newtonian values of T/|W| at neutral points computed by Cutler 
& Lindblom (1992). Two results agree within the maximal errors estimated here except for m = 2 modes 
which are not available for the post-Newtonian calculations. 



Here we have shown that the RCA works rather well in evaluating eigenfrequencies of rotating stars 
in general relativity despite its naiveness. As eigenmode analysis in the RCA is much easier than in 
the post-Newtonian treatment, this approximation can be extensively used in determining oscillatory 
frequencies of realistic neutron stars and/or determining their neutral points. 



T 

Maximal Error of -—— - — Aao ■ 




(1) 



da 



5. Concluding Remarks 



The main drawback of this approximation is that it is hard to include the effect of gravitational radiation 
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which is responsible for damping or growing of corresponding modes. The post-Newtonian treatment can, 
in principle, systematically include this effect by increasing further the order of approximation. In contrast 
the RCA in itself cannot include its effect consistently. Some modification such as an estimation of effective 
damping due to gravitational radiation by using eigenfunctions obtained by the RCA may be needed. 

We are grateful to Dr. Shijun Yoshida for discussions. We would also like to thank the anonymous 
referee for letting us know the investigation by Stergioulas and Friedman, and Prof. John Friedman for 
sending us a copy of their preprint. 

This paper is presented as a part of the Ph.D. thesis of the first author to the Department of Astronomy, 
University of Tokyo. This research has been supported by Grant-in- Aid for Scientific Research of the 
Japanese Ministry of Education, Science and Culture. 
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A. Basic Equations of the RCA 

A.l. Adiabatic perturbations to stationary axisymmetric configurations 

Background spacetime is expressed by the following metric which is used for stationary axisymmetric 
spacetime: 

ds 2 = -e 2v ^dt 2 + e 2a ^ e \dr 2 + r 2 d6 2 ) + e W) r 2 sin 2 6(d<p - w(r, 9)dtf . (Al) 
As we adopt the Cowling approximation here, perturbations to this metric are totally omitted. 
As for the matter variables, perfect fluid is assumed: 

Tab = (p + e)u a u b + pg ab , (A2) 

where u a , e, and p are the 4-velocity, the energy density and the pressure of the fluid, respectively. The 
matter in the equilibrium state is assumed to be polytrope, i.e. p = Ke 1+1 / N where if is a constant. 
Perturbations for the matter are assumed to be adiabatic so that the following Lagrangian relation holds 
for the energy density perturbation Ae and the pressure perturbation Ap: 

— = 7 — > (A3) 
P e 

where 7 denotes the adiabatic exponent. Though this exponent needs not, in general, coincide with the 
equilibrium polytropic exponent, we assume in this paper that they coincide with each other for simplicity. 
Under this assumption, this Lagrangian relation is reduced to the Eulerian one for the Euler perturbations 
for the energy density Se and the pressure dp as follows: 

Sp 5e 

— =7— ■ A4 

P e 

Since the background spacetime has timelike and spacelike killing vectors Jj and all coefficients 
which appear in the perturbed equations do not depend on neither t nor ip. Consequently solutions of the 
perturbed equations have dependence of ~ exp(— iat + imp) on t and ip variables where m is an azimuthal 
eigenvalue and integer. 



A. 2. Introduction of dimensionless variables 

We introduce the following dimensionless variables by using e c and n (variables underlined are 
dimensionless) : 

e = e/e c , (A5) 

p = p/e c = Kelt N -e 1+1 ' N = ne 1+1 ' N , (A6) 

w 4 = w7(w°V^) , (A7) 

a = a/il c , (A8) 

n = n/n c , (A9) 

r = r/r c , (A10) 
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where f2 is the angular velocity. Parameters r c and Q c are defined by: 



ft c = % /4¥i; , (All) 
r c = ^/n c . (A12) 

The quantity ^/k is a dimensionless sound velocity at the center (normalized by speed of light which is 
now chosen as c = 1). This parameter measures the system's strength of gravity and corresponds to the 
expansion parameter in weak gravity limit. 

Further we introduce the Emden function ip, for it will help us to impose the numerical boundary 
condition at the stellar surface properly: 

e = yj N , p = K ^ N+1 . (A13) 



A. 3. Surface-fitted coordinates 

Since a rapidly rotating star deforms from spherical configuration, grid points do not, in general, fall on 
the surface in the ordinary polar coordinates. This makes it difficult to impose the boundary condition at 
the stellar surface. To avoid this difficulty, we introduce the surface-fitted coordinates as in the Newtonian 



analysis (see e.g. Yoshida fc Eriguchi 1995| ). For a configuration with the stellar surface r = R s {9) in the 



polar coordinates (r, 9), new coordinates (f, 6) are defined by 

r = r/R s (0) , = . (A14) 

Although owing to this coordinate transformation the basic equations become slightly complicated, the 
boundary condition at the surface, i.e. vanishing of the Lagrangian perturbation of the pressure 

Ap = 0, (A15) 

can be imposed easily. 

Hereafter, for simplicity, 'tilde' and underlines will be omitted in the equations expressed by this new 
coordinate. 

A. 4. Hydrodynamic Equations in the Cowling approximation 

Basic equations of relativistic hydrodynamics are 

V • (nu) = , (baryon number conservation) (A16) 

VfeT ab = , (energy-momentum conservation) (A17) 
e + p ~ 

de = dn , (the first law of thermodynamics) (A18) 

where n is the baryon number density, and differential operators V and d are the covariant derivative and 
the exterior derivative, respectively. 
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By introducing 3- velocity components on orthogonal basis, 

U = -ie a - v 8v\ 
V = -ire a - u 8v , 



(A19) 
(A20) 
(A21) 



and rewriting the variable Sip as ^ = Sip, the basic equations for the perturbed quantities are written as 
follows. 



For the baryon number conservation, we have, 



(<7 - mfl)$ + 



-0(1 + Kip) v _ 

' NR S (6) C 

— ■0(1 + Kip) - 



d r U 



e v - a d r ip ip(l + Kip)e"- a d r ip 



Rs(6) 



NR S (6) 
ip(l + Kip)e v ~ a 



2d r v + d r a + d r (3 + 



2NAR S (9) 



U 



NR s {6)r 



d B V 



ip(l + KiP)e u ~ a S 
NR s (6)r 



d r V 



+ 



+ 



+ 



e v - a {d e ip-rSd r ip) t ip(l + KiP)e u - a (d e A - rSd r A) 

R s (9)r 
e"- a ip(l + Kip) 



2NAR s {0)r 
((dg - rSd r ){2u + a + (i) + cot I 



NR s {6)r 

0(1 + Kip){u - n)e v+ PR s (6)rs,va.6 
NA 



V 



(a — mQ) 



-0(1 + Kip)e v 



NR s (9)r sin 8 



W = 0. (A22) 



Rclativistic version of the Euler equation is obtained by applying the projection tensor P a b = u a Ub + g a b 
onto the equation of energy-momentum conservation as follows: 



(a-mfl)f/ + [rd r (j + 2((j-n)rd r /3 + 2(uj-n) + 2K(N + l)Brd r ip{uj ~ tt)]e f) ~ a sin6W 

(N + i)ABe-°-f>^ + k(N + l)AB"e-^^-^ = 0, 

(a-mQ)V + [{d g -rSd r )uj + 2{uj-Q)(dg-rSd r )/3 + 2(LU-n)cote 
+2k(N + l)B(dg - rSd r )ip(uj - fi)] e ' 01 sin0W 



- (N + 1)AB- 



-(dg - rSd r )^ + k(N + 1)AB 2 



-(de-rSd r )V = 0, 



(a-mn)W - 



R s (6)r y " K ' R s (0)r 

2(w - tydris - k(lu - nfe^-^Rsiefr 2 sin 2 6 - d r Lu - 2(w - Q) [ d r /3 + - 



(A23) 

(A24) 

e fj - a rsm6U 



- [2(w - n)(dg - rSd r )v - k{(J - n) 2 e 2f3 - 2v R s {8fr 2 sin 2 6{dg - rSd r )uj 

-{dg - rSd r )uj - 2(w - Sl)(d - rSd r f3 + cot 6)} e f3 ~ a sinOV 
+ k<j(N + l)AB(a -{l)e?- 3v R,(6)r8m6* 



+ m(N + l)AB 



- kui{u) - n)e p - 3u R s (6)rsm0 



* = 



R s (6)rsm6 
where functions A, B and S are defined as, 

A = g ao + 2g Qj v j + g jk v j v k = -e 2v + e 2/3 r 2 sin 2 6{u - tl) 2 , 

1 



(A25) 



B = 



1 + Kip ' 



(A26) 
(A27) 
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S{0) 



dR s (9) 



R s (6) dO 



(A28) 
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Fig. 1. — A schematic behaviour of an cigcnfrcqucncy near the neutral point. 'Full calculation' means values 
obtained from full inclusion of gravitational perturbations. 'Correction with maximal error' line is the line 
with the maximal error denoted by vertical arrows and extrapolated from a spherical model. The open circle 
and square denote approximated neutral points. 

Fig. 2. — Critical values of T/\ W\ at neutral points for several orders of f-modes for N = 1 polytrope. In this 
figure, as a measure of strength of gravity the ratio of the mass to the radius in Schwarzschild coordinate 
M/R evaluated for spherical models is adopted. Error bars with open squares for both sides show maximal 
errors expected for the RCA. Integers m are the orders of the mode. Crosses are values of post-Newtonian 
results. Horizontal lines in the upper region denote critical values where equilibrium sequences terminate 
due to mass-shedding. 



Table 1. T/\W\ Values at Neutral Points 



N = 0.5 



N = 1.0 



N = 1.5 



0.001 
0.01 
0.1 
0.5 



m = 2 m — 3 m — 4 m — 2 m — 3 m — 4 

1.36.E-1 9.06B-2 ••• 9.31E - 2 7.90E - 2 

1.28.E-1 8.30B-2 ••• 9.02E - 2 7.39E - 2 

1.47B-1 9.25B-2 6.17E - 2 9.70E - 2 6.95B - 2 4.91B - 2 

7.92B-2 5.11B- 2 3.75B - 2 7.53E - 2 A.WE-1 2.93E - 2 



m — 2 m — 3 m — 4 

5.69B- 2 4.92B-2 

5.66.E- 2 4.84B-2 

5.65B- 2 4.77.E-2 3A4E - 2 

4.14B- 2 2.92B-2 1.98B - 2 



Schematic behaviour of an eigenfrequency near a neutral point 
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